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l^^ ^ Abstract. We establish a connection between finite fields and finite dynam- 

^sj , ical systems. We show how this connection can be used to shed light on some 

problems in finite dynamical systems and in particular, in linear systems. 
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1. Introduction 



00 

q 

" •p^ I There is a natural correspondence between the set GF{p'^) and the set Zp^ of 

Cd ■ r-tuples over Zp, p prime. Furthermore, GF{p'^) is a vector space over GF{p) = 

Zip and hnear transformations over Zp^ correspond to Hnearized polynomials over 

GF(jf). In this ongoing work, we use these facts to study some problems in finite 

dynamical systems and linear automata. 

^ ' In Section 2, we study finite dynamical systems and how our approach can be 

■"ijlj- \ used in the classification problem. In Section 3, we study linearized polynomials 

CO ' and how they can be applied to a problem in linear finite state machines, which in 

"^ \ turn arises from a problem in crystallographic FFTs. 

■^ ' 

o 

^^ \ Finite dynamical systems are important in applications to computational molec- 

r^ ' ular biology. They are useful in microarrays of genes in order to find the best 

"jrt ! model that fits a given data, the so called "reverse engineering problem." (See 



2. Finite dynamical systems 



http://industry.ebi. ac.uk~brazama/Genenetsf) Laubenbacher and Pareigis [2] de- 



fine a finite dynamical system as a function / : fc" ~> /c", constructed by the 
K^ ■ following data: 

k> : 1- fc-={0,f} 

5_j ■ 2. a finite graph F on n vertices. 

C^ I 3. a family of "local" update functions fa : fc" -^ fc", one for each vertex a d F, 

which changes only the coordinate corresponding to a, and computes the binary 
state of vertex a. These functions are local in the sense that they only depend on 
those variables which are connected to a £ F. 

4. an "update schedule" tt, which specifies an order on the vertices of F, repre- 
sented by a permutation tt G 5„. 

The function / is then constructed by composing the local functions according 
to the update schedule vr, that is 

/ = fn{n) O ■ • • O /^(i) : fc" -^ k". 
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In [3], Laubenbachcr and Pareigis called the above function a permutation se- 
quential dynamical system, and extended the definition in different directions. In 
particular, they take k as an arbitrary set. For our purposes it is convenient to 
consider k — Z„i == {0, 1, • • • , m. — 1}, 

Definition 2.1. A finite dynamical system (FDS) is a pair (V, /) where V is the set 
of vectors over a finite field and f : V ^ V. 

Definition 2.2. The state diagram of a FDS (V, /) is the digraph whose vertices are 
members of V and whose edges are the set of all {x, f(x)), where x e V. 

Remark: Note how for an FDS (Zp^ , /), the function can be viewed naturally as an 
r-tuple (/i, /2, • • • , /r) of functions fi : Tip ^ Zp, i < i < r. It is also important to 
note that, using Langrange interpolation, any function from a finite field to itself 
can be realized as a polynomial [4]. Hence each of the fi can be regarded as a 
polynomial over Zp. A similar remark applies to a FDS of the form (Gi^(p"), /). 

Definition 2.3. Two FDSs are isomorphic if their state diagrams are isomorphic. 

Two FDSs are isomorphic if and only if their state diagram are isomorphic as 
digraphs. A limit cycle is simply a directed cycle in the state diagram Sf. A loop 
is a limit cycle consisting of a single vertex and in the case that it occurs, it is a 
fixed point of the FDS /. We denote by Cf the subdigraph of Sf induced by all the 
arcs of the limit cycles. 

Definition 2.4. Let f : V ^ V he a FDS with state diagram Sf and with 
subdigraph Cf of limit cycles. Then x G fc is a vertex in Cf if and only if there 
exists a positive integer m such that /'"(x) = x. The minimum m such that 
/"'(x) = X for all X G Cf is called the order of the system /, denoted by Order(/). 
(See [2]) 

Directed paths in Sf correspond to iterations of / on the element at the beginning 
of the path. Since the set Zp^ is finite, any directed path must eventually enter 
a limit cycle. Thus each connected component of Sf consist of one limit cycle, 
together with transients, that is, directed paths having no repeated vertices and 
ending in a vertex that is part of limit cycle. 

One of the main problems in FDSs is their classification. Loosely speaking, this 
is the problem of determining of two arbitrarily given FDSs whether or not they 
are isomorphic. One of our goals in this work is to facilitate the solution of the 
classification problem though the association given in 

Theorem 2.5. For any fixed basis ai, • • • , a^ of GF{p'^) there is a natural one-one 
correspondence between the FDSs over GF{p^) and those over Zp*". 

Proof. There is a natural correspondence between the sets Zp^ and GF{p^), namely, 
(xi, X2, • • • , Xr) ^-^ xiai + ■ ■ ■ XrUr- Now given / : Zp^ —>■ Zp , define L : GF{p^) -^ 
GF(p'^) such that for each (xi, X2, • • • , Xj.) in Zp^ , L{x) = /i(xi)ai -!-••■ fr{xr)cxr, 
where x = xiai + ■ • • + Xra-r and / — (/i, • • • , fr). Conversely, given L : GF{p^) — > 
GF(jf), if L{xiai + •••Xj.ar) — J/icn + ■■■'UrO^r, there corresponds a function 
/ == (.fi, ./2, • • • , fr) such that /(xi, • • • ,Xr) = (j/i, • • • , 2/r), where fi{xi) = yt for 
each i — 1,2, ■ ■ ■ r. Since this correspondence is onto and the two sets are finite with 
the same number of elements, it is also one-one. 
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Corollary 2.6. If 5*1 = {Zp^ , f) and 5*2 = (Zp^ , f) are FDSs and / corresponds to 
L{x) with respect to the basis ai, • • • ,ar and /' corresponds also to L{x), but with 
respect to another basis, then Si and 5*2 are isomorphic. 

This latter corollary says that our approach is quite useful for the classification 
problem. On the other hand, 

Theorem 2.7. For any fixed basis ai, ■ ■■ , a^ of GF(p^), there is a natural corre- 
spondence between the FDSs over {GF(p^))'^ and those over (Zp)*"". 

In other words, it is redundant to study both types of these FDSs, but each is 
important given the classification problem. 

3. Linear FDSs and linearized polynomials. 

A linear finite dynamical system or a linear (autonomous) finite state machine 
is a FDS (Zp^ , /) in which / is a linear transformation on Zp^ regarded as a vector 
space over Zp. We shall see in this case that there is a useful correspondence between 
linear FDSs and linearized polynomials. 

The correspondence x ^ x^ , z = 0, 1, • • • ,p^~^, gives the Galois automorphisms 
of GF{p^). A linearized polynomial L{x) is a polynomial generated by these auto- 
morphisms. In other words, L{x) = J2l=o ^t^^ j where Ai £ GF{p^.) We note that 
iiy,z€ GFip"-) and A e GF{p), then L{x + y) = L{x) + L{y) and L\x) = \L{x). 
Thus, L{x) is a linear function on GFijf) regarded as a vector space over GF{p). 
Furthermore the correspondence between GF{p^) and Zp^ given in Theorem 1 is 
an isomorphism as a vector space over Zp . Since there are {p''Y linearized poly- 
nomials, this coincides with all the linear functions on Zp^ . It is easy to sec that if 
/ : Zp' -^ Zp^ is a FDS associated to the linearized polynomial L{x) = X][=o ^i^^ i 
then kcr/ is the set of all roots of L{x). So, / is invertible if and only if the only 
root of L in GF(p'^) is 0. 

Given a linear autonomous machines S = {Zp^,F), if / is a nonsingular linear 
transformation, i.e., an invertible matrix over Zp, then the state space Zp*^ decom- 
poses into disjoint "orbits" or "cycles." Based on the above observation, the same 
also holds for machines {GF{p^), L) and it is interesting to note how properties of 
linearized polynomials determine this orbit structure, much in the same way as the 
properties of the elementary divisors of / determine the orbit structure of (Z^, /) in 
the classical theory. However, our motivation for studying linearized polynomials 
stems from a more general problem which arises in crystallographic FFTs [5]. Let 
us briefiy describe this problem. 

Crystallographic data can introduce structured symmetries in the inputs of a 
multidimensional discrete Fourier transform, which in term introduce symmetries 
into the outputs. In order to avoid redundant calculations, it is of interest to 
exploit these symmetries. Assuming that symmetries are given by an n x n matrix 
S over Zp, for prime p edge length, we can reduce the complexity of the FFT by 
determining a matrix M with MS = SM and M*5 =: 5'M* (where M* denotes 
the transpose of M) that minimizes the number of "MS'-orbits." A vector x £ Z" 
belongs to an MS-orbit of length k if and only if M'^x = S'^x for some i. The cases 
n — 2 and n — 3 are of particular interest. 

For n = 2, for example, i.e., F : Zp x Zp ^ Zp x Zp, this corresponds to the 
study of linearized polynomials F{x) = Ax'p + Bx, where A,Bg GF{p^), where F 
is an invertible map. The first question is, when is F{x) invertible in GF{p^)l 
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Lemma 3.1. F{x) is invertiblc if Ap+^ 7^ BP+\ 

Proof. F{x) is invertiblc if and only if it is one-one, i.e., if ker _F = 0. In other 
words, if F{x) ~ has only a; = as a solution over GF(p^). But x ^ and 
AxP + Bx = imply that X^^^ = — ^ and so raising both sides to the power p+ 1, 

we obtain xp'-^ == |^ and Ap+^ = BP+\ 

In the remainder of this section we consider the class Cp of linearized polynomials 
L{x) = J2l=o AixP where Ai e GF{p). These types of polynomials have important 
properties which we outline below. 

Property I. If L{x),L'{x) G Cp then L{L'{x)) = L'{L(x)). (Note that this means 
that the corresponding matrices commute). 

Property II. Given L(x) G Cp, the class of L'{x) satisfying Property I is precisely 
r 

Property III. Using a normal basis, the matrix corresponding to L{x) G Cp is 
symmetric. Consequently, the transpose matrix also commutes. 

Definition 3.2. If L{x) = X]I=o ^j^^ i^ ^ linearized polynomial, then its associate 

isl{x)^Y.ZlA^X\ 

Definition 3.3. U{x) = L{x) and L"+i(a;) = L(L"(a;)). That is, L"(a;) is the n-fold 
composition of L with itself. 

Property IV. L^{x) = x modulo xP = x if and only if (l{x))"' = 1 modulo x'' = 1. 

4. Systems over Zp,^ 

In this section we study FDS over Zp^, that is systems (V, /) where / : Zp^^ -^ 
ZprJ and V is the Zp^-module Zp^''. We use Definitions 2.2, 2.3, and 2.4, but now 
we are working in the Zp^-module V . Therefore a linear FDS is an endomorphism 
of the Zpi. -module V . 

Let g : Zp" -^ Zp^ be a bijection. Then the product function g^ : Zp™' -^ Zp„ 
given by g{ai,a2,- ■■ ,ar) ^ (g(ai), 5(02), • • • ,3(0^)) is a bijection too. 

Proposition 4.1. Let / : Zp.'' ^ Zp.'' be a FDS. Let / : Zp"'' -^ Zp"'' be the FDS 
such that g^ o f = f o g^ ■ Then / and / have state diagrams isomorphic. 

Proof. Since g^ is a bijection, there exists {g"^)^^ = g^^ . Then the system /i = 
g^'' o f o g'' has the same state diagram to /. In fact, set xi ~ g^^{x) and yi = 
9~^{y) — (5 '^ ° / ° 9^){xi)- Now, suppose (x, y = f{x)) is an edge in the state 
diagram of /. Then {xi,yi) is an edge in the state diagram of /i. On the other 
hand, g^ ° fi = f ° g^ ■ So, f — fi and our claim holds. 

Definition 4.2. With the notation above, if / is a linear FDS over Zp^ then the 

system / will be called the linear system associated to / by the bijection g. 

We use the the linear system associated to a non-linear FDS / to describe its state 

diagram and its order. 

Example 

Let / : (Z23)^ -^ (^23)^ be a linear system, given by f{a,b) — (5b, a + 26). For 
any bijection g : (Z2)'^ -^ Z23 we have an induced systems / with the same state 
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diagram of the system /. Now using the method given in [1], we find the Order of 
/ and the Order of / for any bijection g. The matrix 

has minimal polynomial 

m{x) = {x - 1)^ 
and the Order of A modulo 2 is e = 2. Since 

5 2 






^ j ( mod 23 ), 



the largest positive integer /3 such that A^ = I { mod 2'') is /3 = 1. Then A^ has 
Order 4 and A has Order 8 modulo 2^. 

5. Future Work 

We will exploit the ideas presented here to seek a polynomial solution to the 
reverse engineering problem. We also make use of our theory to develop an efficient 
algorithm to determine, given a matrix S of symmetries, a matrix M that minimizes 
the number of MS'-orbits in the precomputation phase of crystallographic FFTs. 
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